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Abstract 

O ' 

' The dynamics of test particles in f(Q) modified Gauss-Bonnet gravity is investigated. It is shown that 

in f(G) gravity models with non-minimal coupling to matter, particles experience an extra force normal to 
their four- velocities and as a result move along non- geodesic world- lines. The explicit form of the extra 
force depends on the function of the Gauss-Bonnet term included in coupling term. The effects of this force 
on the relative accelerations of particles are studied. 
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1 Introduction 

J> ■ The f(G) modified Gauss-Bonnet gravity (or f(Q) gravity for short) is among the plethora of models proposed 

recently to explain certain physical observations such as the late time accelerated expansion of the universe [1] . 
It is well known that adding the Gauss-Bonnet curvature Q to the usual scalar curvature in Einstein-Hilbert 
action does not alter the equations of motion of general relativity due to the fact that it is a topological 
invariant in four dimensions, even though it may have some other non-trivial contributions [5J [3J. However this 
term can be incorporated into space-time dynamics in several ways, namely by considering higher dimensions, 
e.g as in [4] or via non-minimal coupling to some scalar field in arbitrary space-time dimensions, which is 
also called modified Gauss-Bonnet gravity, as in string gravity, see e.g [6]. In f{Q) gravity an arbitrary 
function of the Gauss-Bonnet curvature is added to the Einstein-Hilbert action resulting in modified equations 
of motion. Physical and astrophysical implications of such a modification have been extensively studied, e.g. 
in [3 [H El [TUl HU HH H3l HH [JSl US] . An extension of the model has been considered in [H [H DH [50] • A 
review and comparison with f(R) gravity model may be found in [5T] . 

In the present work we aim to investigate the dynamics of test particles in f(G) gravity models. This is 
motivated by the essential role of particle dynamics in physics and astronomy. In particular, in the present 
context, the particle dynamics may be used to put constraints on the form of the Lagrangian as it has been 
done in [55] for the case of R n gravity. We show that within the usual f(Q) gravity, i.e. f(G) gravity minimally 
coupled to matter, test particles move along the background geodesies, but if we consider non-minimal coupling 
with matter extra force arises due to the non-minimal coupling resulting in non-geodesic motions. As it has 
been shown in [23) , a similar effect can emerge in models in which an arbitrary function of the scalar curvature 
is non- minimally coupled to matter in the framework of f (R) gravity [M] [55] . It should be noted that a non- 
minimal coupling does not necessarily result in non- geodesic motion. This point has been shown in [26j for a 
i?-matter coupling with various choices of the matter field. 

In this work we consider an action which consists of the usual ,f{Q) gravity action and an arbitrary function 
of Q non- minimally coupled to matter Lagrangian. This is a slight modification to the model proposed in 
[571 [55] which has attracted some attention, namely in [55] in which a modified Gauss-Bonnet model with 
non-minimal coupling has been used to investigate the cosmic acceleration and crossing of phantom divide. 
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In what follows we first give a brief review of the f(G) gravity and by deploying the energy-momentum 
tensor of a perfect fluid we show explicitly that within the framework of f(G) gravity test particles move along 
the background geodesies. In the next section we show that in a f(G) gravity with non- minimal coupling 
to matter an extra force will be exerted on particles resulting in non- geodesic motion. Then we obtain the 
world-line deviations of particles as a generalization of the geodesic deviation equation for both f(G) and f{R) 
gravities with non-minimal coupling and compare them. We conclude the work with a discussion of the results. 



2 Motion in f(Q) gravity 

In f(G) gravity, the space-time equation of motion can be obtained from the action [1] 

S = jV^g QpR + f(G) + L^j d A x (1) 

where k 2 = 8irG is the Einstein-Hilbert coupling constant which we set to unity henceforth, L is the matter 
Lagrangian, and f(G) is an arbitrary function of G = R 2 - 4i? Q/3 i? Q/3 + R a p K x R°P kX . Varying this action with 
respect to metric components g^ v leads to the equation of motion 

T» v = G» v - <T7 + Af'H^ (2) 

where G» v is the Einstein tensor, T» v = -L- g 5{L S ^ 9) is the energy-momentum tensor, / and /' stand for f(G) 
and J ^ respectively, and 

+2R^ al3 "R aP - 2G t " y V 2 - flVV 
-2.g^i? a/3 V Q V /3 + 2R av V Q V^ 

+2R" a V a V u -2R" al3u V a V ( 3. (3) 

Now we show explicitly that within the framework of f(G) gravity test particles move along the background 
geodesies. Taking the covariant derivative of both sides of equation ([2]) we can simplify the result to 

V„T^ = - V/ (4) 

where the right hand side vanishes identically. For a perfect fluid with pressure p and energy density p we have 

= P h^ + pu"u u (5) 

where W = + u^u v . Inserting this into ^ and projecting the resulting equation parallel and normal to 
results in the continuity equation 

P = -(p + p)e (6) 

where 9 = W v u v , and the equation of motion 

Vp + u^p +{p + p)^ = (7) 

where an over-dot means ii Q V a and the normalization u M it M = — 1 has been used. For a perfect dust we have 
p — 0, p — nm where n is the particles number density and m is the mass of each particle and thus the last 
equation reduces to 

U» = (8) 

corresponding to a geodesic trajectory. For massless particles we have u^u^ = and the above procedure again 
results in geodesic equation, as in general relativity. 
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3 Non-geodesic motion 

In this section we consider an extension of the above model in which an arbitrary function of the Gauss-Bonnet 
curvature is non-minimally coupled to matter. The action may be written in the following form 



S = J^g QiZ + f(G) + (1 + XF(G))Lj d A 2 



(9) 



in which F(Q) is an arbitrary function of Q and A is a coupling constant. This may be compared with the 
action 

^R + fiG^^x (10) 

introduced in [27J HE] • The equation of motion resulting from variation of the action © reads 

(1 + \F)T» V = - g^f + AH^(f + XLF') (11) 
where F = F(Q) and F' = dF d g ■ Taking the divergence of both side of this relation we obtain 

V„T^ = -2^—(( g ^L - T» V )V V G + K^VyL) (12) 

1 + AT 

where K ^ = &RW V + AR^ aPl R ua ^ - WR£R av + 8i?' ia/3l/ i? Q/3 . The above equation can be used to obtain the 
matter equation of motion. For the perfect fluid described by energy-momentum tensor given in ([5]) projection 
parallel and normal to results in 

A F' / \ 
P+(p + P )9 = — ((L|p)(j|FV^) (13) 

and 

XF'W 

(p + pW + v»p + u»p = —^L.((L-p)^g 

+K l/a V a L) (14) 

respectively. In the absence of non-minimal coupling the right hand side of this equation vanishes and we 
arrive at the geodesic equation. The right hand side of equation (|14j) is non-vanishing in general which means 
that test particles feel an extra force which depends on F(Q). This extra force which originates from the 
non-minimal coupling is normal to the trajectories of test particles and results in non-geodesic motions. 
For a perfect dust with L = —p, equation (fi"4"f reduces to 

tf = Y^(y v g + K va v a hL P ). (15) 

For the low energy effective action with matter part given by 

L=p(4>,X) (16) 

in which is a scalar field and 

x = iv^vv 

the energy momentum tensor may be transformed into the perfect fluid form given by ^ by setting 



provided V M </> being time-like [5U]. The equation of motion can then be obtained from (TT4"]) . 
Equation (|12[) may be compared with the one obtained in [23j for f(R) gravity 

V V T" V = A/ f (g» v L - T^ V )V„R. (17) 

1 + A/2 

In the latter the perfect fluid Lagrangian L = p results in a vanishing extra force but equation (|14p shows that 
this is not the case for f(Q) gravity. Similarly by choosing L = —p, the right hand side of ([13]) does not vanish 
but it vanishes for the counterpart equation in f(R) gravity. 
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4 World-line deviations 



In general relativity the equation of geodesic deviation is used to describe relative motions of test particles 
in a frame independent way. Here we obtain its counterpart in f(G) gravity with non- minimal coupling to 
matter which we call the equation of world-line deviation. To this end, we consider a congruence of world-lines 
parametrized by e and compare the equations of motion (|15p for two particles moving on nearby trajectories 



x^(t) and ie m (t) + en M (r) respectively, where e is small and n M (r) = dx ' Jf^ . On keeping only terms which are 
linear in e at most and defining 

— \F' 

we obtain 

n» = -R» av /3U a n u u + h£n a V a F' + JvKn" + u"n v ) (18) 



where we have used 

Du^ Dn^ 



De Dt 



(19) 



= + R" w pu a n v uf l . (20) 

De Dt Dt De 

The first terms in the right hand side of equation (fT5)) shows the usual effect of the background curvature as 
in the well-known geodesic deviation equation. The second term shows a contribution from the non-minimal 
coupling which is normal to the four- velocity and depends on n M . The last term is proportional to the relative 
velocity. Such velocity-dependent terms also appear in the Palatini formulation of the f(R) extended gravity 
without coupling to matter |31j . The spacial components of (| 18|) may written in the form of 



a = ai +a t (21) 

where l,t stand for longitudinal and transverse with respect to the velocity. The transverse component is 
responsible for the deviation from the geodesic path. 

By following the procedure described above, one can obtain a similar deviation equation for f(R) gravity 
with non-minimal coupling. Starting from equation (|17p we can obtain the following equation of motion for 
perfect dust particles 



h»V v R (22) 



l + A/ 2 

from which we arrive at equation (|18p again, but with T v replaced by 

-Afo 

An interesting difference between T% and T v is that aside from a coupling constant A the former depends only 
on purely geometric objects but the latter depends on both geometric objects and matter mass density directly. 
By an expansion of the form F(Q) = F(0) + F'(0)Q + ■ ■ •, we obtain 

T v = -\F'(0)(V v g + K va V a lnp) + • • • 

which is quadratic in curvature. A similar expansion for the case of f(R) gravity with non-minimal coupling 
results in an expression which is linear in curvature. From this vantage point, the deviation from geodesies is 
weaker in the model considered here. For large curvatures the deviation would be stronger in non- minimal 
/((?) gravity. 

Equation (|18|) may be used to study the stability of orbits of particles in the framework of the model under 
consideration. Given a specific orbit of a particle in a space-time which is a solution of the equation of motion 
(fTTjl (or its f(R) gravity counterpart), one can solve equation (fT8| for the deviation from the fiducial orbit. 
This can then be put into 

A=limiln(vJ22^ (23) 



r — >oc f 
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to compute the relevant Lyapunov index which in turn is used to determine the stability behaviour of the orbits 
(see e.g. |32j for a discussion in the context of general relativity). Non- minimal coupling contributes to this 
through ^''-dependent terms. 

5 Discussion 

In this work we showed that in f(G) gravity model with non- minimal coupling of an arbitrary function of the 
Gauss-Bonnet curvature to the matter Lagrangian, test particles move along non-geodesic trajectories due to 
the presence of extra forces originated from the coupling. These extra forces are perpendicular to the particles 
velocities and depend on both the choice of the function couples to matter and the matter Lagrangian. These 
results qualitatively agrees with those obtained in |23j for the case of non- minimal coupling in f(R) gravity. 
However, as we have shown, there are important distinctions between these models, namely the extra force in 
non-minimal f(R) model vanishes for certain usual choices of matter Lagrangian but is does not in non-minimal 
f(G) model. We also generalized the geodesic deviation equation to a world-line deviation equation for both 
f(G) and f(R) gravities with non-minimal couplings. According to this world-line deviation equation, relative 
accelerations of test particles depends on background curvature (as in general relativity), relative distances, 
and relative velocities of the particles. This is similar to the results obtained in the case of f(R) gravity in 
its Palatini formulation without non-minimal coupling even though their origins are different. We showed 
that compared with non-minimal f{R) gravity the deviations from geodesies is weaker in non-minimal /(<?) 
gravity for small curvatures. It was also shown schematically how the world-line deviation equation we have 
obtained might be deployed to study the stability of orbits of particles in the extended gravity model considered 
here. More elaborate study of these relative motions would be possible by taking a generalized Raychaudhuri 
equation into account. The issue of energy conditions and stability would also be of interest in the context of 
the model considered here. 
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